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ABSTRACT 

We explore the implications for the high redshift universe of "state-of-the-art" 
models for the chemical and spectrophotometric evolution of spiral galaxies. The mod- 
els are based on simple "scaling relations" for discs, obtained in the framework of Cold 
Dark Matter models for galaxy formation, and were "calibrated" as to reproduce the 
properties of the Milky Way and of nearby discs (at redshift z ~ 0). In this paper, we 
compare the predictions of our "hybrid" approach to galaxy evolution to observations 
at moderate and high redshift. We find that the models are in fairly good agreement 
with observations up to z while some problems appear at higher redshift (provided 
there is no selection bias in the data); these discrepancies may suggest that galaxy 
mergers (not considered in this work) played a non negligible role at z >1. We also 
predict the existence of a "universal" correlation between abundance gradients and 
disc scalelengths, independent of redshift. 
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1 INTRODUCTION 

Important progress has been made in the past few years in 
our understanding of galaxy evolution, mainly due to the 
number and the quality of observations of galaxies at in- 
termediate redshifts (up to z ~1). Observations of ground- 
based telescopes, combined with results from the Hubble 
Space Telescope, provided photometric and kinematic data 
on the size, luminosity and rotational velocity of galactic 
discs as a function of redshift (Lilly et al. 1998, Schade et 
al. 1996, Vogt et al. 1997, Roche et al. 1998, Simmard et al. 
1999). Although the question of various selection effects is 
not settled yet, these data allow, in principle, to probe the 
evolution of normal galaxies and to test theories of galaxy 
formation. 

Several attempts have been made to interpret these data 
in the framework of currently popular models of hierarchical 
growth of structure in the Universe. Combined with simple 
assumptions about angular momentum conservation, such 
models lead to simple scaling relations as a function of red- 
shift and may successfully reproduce some of the aforemen- 
tionned high redshift data (Mo et al. 1998, van den Bosch 
1998, Steinmetz and Navarro 1999, Firmani and Avila- Reese 
2000, Avila-Reese and Firmani 2000). This "forwards" ap- 
proach infers the structural properties of discs from those 
of the corresponding dark matter haloes and does not con- 
sider in detail the driving force of galaxy evolution, namely 
the transformation of gas to stars; for that reason, assump- 
tions have to be made about the variation of the mass to 



light ratio and the colours of a disc as function of its mass 
and of redshift, while the associated chemical evolution is, 
in general, ignored. 

On the other hand, several studies (Ferrini et al. 1994, 
Prantzos and Aubert 1995, Chiappini et al. 1997, Boissier 
and Prantzos 1999) focused mainly on the properties of lo- 
cal discs, for which a large body of observational data is 
available. In those studies, the observed chemical and pho- 
tometric profiles of discs are used to constrain the radial 
variation of the Star Formation Rate (SFR) and of the infall 
rate. These multi-zone models naturally lead to a relation 
between luminosity and size (i.e. exponential scalelength) of 
discs. The fact that the Milky Way is a typical disc galaxy 
helps to "calibrate" such models, but there is no simple pre- 
scription as to how to extend them to other discs. 

In principle, one can utilise a "backwards" approach 
and try to infer the properties of high redshift discs from 
the histories of the models that reproduce well the local disc 
population. Such an approach has been adopted in Cayon et 
al. (1996), Bouwens et al. (1997), Prantzos and Silk (1998) 
and Bouwens and Silk (2000). In the latter work, a compar- 
ison is made between simplified versions of the "forward" 
and "backwards" approaches. 

A kind of "hybrid" approach is adopted in Jimenez et 
al. (1998): they relate the disc surface density profile to the 
properties of the associated dark matter halo, while they 
adopt radially varying SFR and infall timescales such as to 
reproduce in detail current profiles of the Milky Way disc. 
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While Jimenez et al. (1998) applied their model to the study 
of Low Surface Brightness galaxies, we developed a modi- 
fied and more detailed version of this "hybrid" model and 
applied it extensively to the study of local discs. In a series 
of papers (Boissier and Prantzos 1999, 2000; Prantzos and 
Boissier 2000; Boissier et al. 2001) we showed that this model 
can readily reproduce a large number of local disc prop- 
erties: Tully-Fisher relations in various wavelength bands, 
colour-colour and colour-magnitude relations, gas fractions 
vs. magnitudes and colours, abundances vs. local and inte- 
grated properties, abundance gradients, as well as integrated 
spectra as a function of a galaxy's rotational velocity. A cru- 
cial ingredient for the success of the model is the assumption 
that the infall rate scales with mass of the galaxy, i.e. mas- 
sive discs form the bulk of their stars earlier than low mass 
ones. 

Stimulated by the success of this "hybrid" model in re- 
producing properties of the local disc population, we explore 
in this work its implications for the high redshift Universe, 
comparing its predictions to data of recent surveys. The plan 
of the paper is as follows: In Section 2, we present the mod- 
els that were developped in the previous papers of this series 
in order to match the properties of the Milky Way and of 
nearby spirals; we also present the adopted probability dis- 
tributions of the two main parameters of our models, namely 
the disc rotational velocity Vc and its spin parameter A. The 
predictions for the evolution of spirals, deduced from those 
models, is shown in Section 3. In Section 4 we compare the 
results of models with observations up to z = 1. In Sec- 
tion 5, we exploire the implications at higher redshifts and 
we compare our results to scalelength distributions observed 
up to z = 3.5; we also predict that abundance gradients and 
scalelengths should present at all redshifts the same corre- 
lation as the one observed locally. The conclusions of our 
"backwards" approach for the xploration of the high red- 
shift Universe of disc galaxies are summarized in Section 6. 



2 MODELS AND RESULTS FOR NEARBY 
DISCS 

In this section we briefly describe the main ingredients and 
the underlying asumptions of our model for the chemical and 
spectro-photometric evolution of spiral galaxies. The model 
has been "calibrated" in order to reproduce the properties 
of the Milky Way disc (Sec. 2.1), then extended to other 
spirals with the help of simple scaling relations that allow 
to reproduce fairly well most of the observed properties of 
nearby spirals (Sec. 2.2). In Sec. 2.3 we present the adopted 
distribution functions concerning the two main parameters 
of our disc models, namely the rotational velocity Vc and 
the spin parameter A. 

2.1 The Milky Way model 

The model for the Milky Way disc is presented in detail in 
Boissier and Prantzos (1999, hereafter Paper I). The galactic 
disc is simulated as an ensemble of concentric, independently 
evolving rings, gradually built up by infall of primordial 
composition. The chemical evolution of each zone is followed 
by solving the appropriate set of integro-differential equa- 
tions (Tinsley 1980), without the Instantaneous Recycling 



Approximation. Stellar yields are from Woosley and Weaver 
(1995) for massive stars and Renzini and Voli (1981) for in- 
termediate mass stars. Fe producing SNIa are included, their 
rate being calculated with the prescription of Matteucci and 
Greggio (1986). The adopted stellar Initial Mass Function 
(IMF) is a multi-slope power-law between 0.1 Mq and 100 
Mq from the work of Kroupa et al. (1993). 

The spectrophotometric evolution is followed in a self- 
consistent way, i.e. with the SFR ^(t) and metallicity Z(t) 
of each zone determined by the chemical evolution, and the 
same IMF. The stellar lifetimes, evolutionary tracks and 
spectra are metallicity dependent; the first two are from the 
Geneva library (Schaller et al. 1992, Charbonnel et al. 1996) 
and the latter from Lejeune et al. (1997). Dust absorption 
is included according to the prescriptions of Guiderdoni et 
al. (1998) and assuming a "sandwich" configuration for the 
stars and dust layers (Calzetti et al. 1994). 

The star formation rate (SFR) is locally given by a 
Schmidt-type law, i.e proportional to some power of the gas 
surface density E 9 and varies with galactocentric radius R 
as: 

V(t,R) = aZ g {t,R) 15 V(R) R' 1 (1) 

where V(R) is the circular velocity at radius R. This radial 
dependence of the SFR is suggested by the theory of star 
formation induced by density waves in spiral galaxies (e.g. 
Wyse and Silk 1989). The efficiency a of the SFR (Eq. 1) is 
fixed by the requirement that the observed local gas fraction 
a g (Ro=8 kpc)~0.2 is reproduced at T=13.5 Gyr. 

The infall rate is assumed to be exponentially decreas- 
ing in time with a characteristic time r. In the solar neigh- 
borhood we adopt t~7 Gyr in order to reproduce the local 
G-dwarf metallicity distribution (Paper I). In order to mimic 
the "inside out" formation of the disc, r is assumed to be 
shorter in the inner zones and larger in the outer ones. 

The really "free" parameters of the model are the ra- 
dial dependence of the infall timescale t(R) and of the SFR 
ty(R)- It turns out that the number of observables explained 
by the model is much larger than the number of free pa- 
rameters. In particular the model reproduces present day 
"global" properties (gas, O/H, SFR, and supernova rates), 
as well as the current disc luminosities in various wavelength 
bands and the corresponding radial profiles of gas, stars, 
SFR and metal abundances; moreover, the adopted inside- 
out star forming scheme leads to a scalelength of ~4 kpc in 
the B-band and ~2.6 kpc in the K-band, in agreement with 
observations (see Paper I). 

2.2 Extension to other disc galaxies 

In order to extend the model to other disc galaxies we adopt 
the "scaling properties" derived by Mo, Mao and White 
(1998, hereafter MMW98) in the framework of the Cold 
Dark Matter (CDM) scenario for galaxy formation. Accord- 
ing to this scenario, primordial density fluctuations give rise 
to haloes of non-baryonic dark matter of mass M, within 
which baryonic gas condenses later and forms discs of maxi- 
mum circular velocity Vc- It turns out that disc profiles can 
be expressed in terms of only two parameters: rotational 
velocity Vc (measuring the mass of the halo and, by assum- 
ing a constant halo/disc mass ratio, also the mass of the 
disc) and spin parameter A (measuring the specific angular 
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momentum of the halo). If all discs are assumed to start 
forming their stars at the same time, the profile of a given 
disc can be expressed in terms of the one of our Galaxy (the 
parameters of which are designated hereafter by index G): 

Ra A V C 

RdG Ag VcG 

and 

So = (A-Y 2 

Sog V Ag / Vcg 

Eqs. 2 and 3 allow to describe the mass profile of a 
galactic disc in terms of the one of our Galaxy and of two 
parameters: Vc and A. The range of observed values for the 
former parameter is 80-360 km/s (with Vcg=220 km/s), 
whereas for the latter numerical simulations give values in 
the 0.01-0.15 range (with Ag ~0. 03-0. 06). Larger values of 
Vc correspond to more massive discs and larger values of A 
to more extended ones. Although it is not clear yet whether 
Vc and A are independent quantities, we treated them as 
such and we constructed a grid of 25 models caracterised by 
80 < Vc < 360 km/s and 1/3 < A/A G < 3. 

As discussed in Boissier and Prantzos (2000, hereafter 
Paper II) the resulting disc radii and central surface bright- 
ness are in excellent agreement with observations, except for 
the smallest values of A (~ 1/3 Ag); this "unphysical" value 
leads to galaxies ressembling to bulges or ellipticals, rather 
than discs. 

The two main ingredients of the model, namely the Star 
Formation Rate ^f(R) and the infall time-scale t(R), are af- 
fected by the adopted scaling of disc properties in the fol- 
lowing way: 

- For the SFR we adopt the prescription of Eq. 1, with 
the same efficiency a as in the case of the Milky Way (i.e. the 
SFR is not a free parameter of the model) . In order to have 
an accurate evaluation of V(R) across the disc, we calculate 
it as the sum of the contributions of the disc and of the dark 
halo, the latter having a density profile of a non-singular 
isothermal sphere (see Paper II). 

-The infall time scale is assumed to decrease with both 
surface density (i.e. the denser inner zones are formed more 
rapidly) and with galaxy's mass, i.e. r[Afa, E(-R)]. In both 
cases it is the larger gravitational potential that induces a 
more rapid infall. The radial dependence of r on is 
calibrated on the Milky Way, while the mass dependence of 
t is adjusted as to reproduce the properties of the galactic 
discs. 

This model, "calibrated" on the Milky Way and hav- 
ing as main parameter the infall dependance on galaxy 
mass, reproduces fairly well a very large number of prop- 
erties of spiral galaxies at low redshift (Paper II): Tully- 
Fisher relations in various wavelength bands, colour-colour 
and colour-magnitude relations, gas fractions vs. magnitudes 
and colours, abundances vs. local and integrated properties, 
as well as spectra for different galactic rotational velocities. 
The main assumption of the model is that infall onto massive 
discs occured earlier than infall onto low mass galaxies. More 
recently, in paper IV (Boissier et al., 2001) we used a homo- 
geneous set of observational data (a subset of the data pre- 
sented in Boselli et al., 2000) in order to test this hypothesis. 
We determined the gas fraction <jq = Ma as /Mtot and the 
star formation efficiency e in a large number of "normal" spi- 
rals as a function of circular velocity Vc- We found that e 



is independent of Vc while ctg decreases with Vc- Taken at 
face value, those findings imply that low mass discs had not 
access to their gas reservoir as early as massive ones, oth- 
erwise they should have already turned most of their gas 
into stars (because their current star formation efficiency is 
similar to the one of massive discs, and we assume that this 
has always been the case). This crucial finding justifies the 
assumptions we made concerning the form of the infall. 

We notice that in another recent work, Brinchmann and 
Ellis (2000) found that the total stellar mass density in mas- 
sive galaxies is constant over the redshift range 0.2 < z < 1. 
This implies that the stellar content of these galaxies must 
have formed at higher redshifts, which is in qualitative agree- 
ment with our assumptions. 

Finally, it is worth noticing that the properties of abun- 
dance and colour gradients among local spirals are readily 
reproduced by the generalization of the Milky Way model 
and the assumed radial dependence of the SFR and infall 
timescales (Prantzos & Boissier 2000, herefter Paper III). 



2.3 Distributions of the parameters Vc and A 

Our models presented in Sec. 2.2 simulate the evolution of 
individual disc galaxies, characterised by Vc and A. In order 
to treat the evolution of the disc population, we need to 
make some assumptions about the probability distributions 
of the relevant parameters. In this paper we assume that the 
distributions Fa (A) and Fv(Vc) are independent and do not 
evolve in time. We adopt the velocity distribution suggested 
in Gonzalez et al. (2000): 



F v (V)dV = *« 



-V 



VJ 



exp 



dV 
V* ' 



(4) 



The parameters V*, (3 and n are determined in Gon- 
zalez et al. (2000) on the basis of observed Tully-Fisher re- 
lationships and luminosity (Schechter-type) functions. We 
adopt here the set of parameters of their Table 4 (fifth row, 
LCRS-Courteau data) corresponding to the velocity inter- 
val covered by our models. By construction, the distribution 
of Eq. (4) is normalised to the local luminosity function for 
spirals. 

The distribution of the spin parameter is given by 



F\(X)d\ = 



2tto\ 



exp 



to 2 (A/A) 



<2o\ 



dX 

T 



(5) 



with A=0.05 and a\=0.5, according to numerical simulations 
(see e.g. MMW98) . The distribution Fx (A) is normalised to 
unity. 

Distributions F\(X) and FyiYc) are shown in Fig. 1. 
Low mass discs are more numerous than high mass ones, 
while the spin parameter distribution favours discs with in- 
termediate A values. Notice that in this work we study discs 
with Vc >80 km/s and A <0.1, since discs with larger A val- 
ues lead to Low Surface Brightness galaxies (e.g. Jimenez et 
al. 1998), which will be the subject of a forthcoming paper 
(Boissier and Prantzos, in preparation). 

In the following, we assume that the number of 
discs with velocity Vc and spin parameter A is dN — 
Fv{Vc)F\(X)dVcdX; we compute then the distribution cj>(Q) 
of various quantities Q (magnitude, radius, etc) , defined by 
dN = (p(Q)dQ. 
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Figure 1. Adopted distributions of the two parameters carac- 
terising our models of disc galaxies: spin parameter A (left) and 
circular velocity Vc (right). The A distribution is normalised to 
J F(\)d\=l, while the Vc distribution is normalised to the local 
luminosity function in Gonzalez et al. (2000). See text (Sec. 2.3) 
for details. 

We make the assumption that the circular velocity and 
spin parameter of a galaxy are determined at its formation 
and do not evolve in time. Since the distributions Fv(V) and 
Fa (A) are always the same, any evolution in <f){Q) results 
from changes in intrinscic properties of individual galaxies, 
according to our models. It is clear that this assumption 
is an (important) oversimplification, since some spirals may 
merge (into e.g. ellipticals), modifying Fv(V). The history 
of spirals that we present concerns only fairly isolated dies, 
that have not suffered major merger episodes. 



3 MODEL PREDICTIONS FOR THE 
EVOLUTION OF SPIRALS 

In Fig. 2 we present the results of our models concerning the 
evolution of the distribution functions of various quantities: 
disc scalelengths Rb in the B-band (top left), central surface 
brightness fio(B) in the B-band (top right), B- magnitude 
Ms (middle left), colour index (U-V)ab (middle right), av- 
erage oxygen abundance in the gas [O /H] (bottom left) and 
star formation rate ip (bottom right). The corresponding dis- 
tributions are shown at various redshifts, starting at 2=4.2 
(dashed curves) and ending at 2=0 (solid curves). A cos- 
mological model with matter density f2M=0.3, f2A=0.7 and 
Ho =65 km/s/Mpc is adopted throughout this work. 

The comparison of our results to observations of local 
spirals has been done in previous papers (Papers II and III) . 
Here we comment on the evolution of the distribution func- 
tions <j>(Q) that we obtain. As can be seen from Fig. 2, there 
has been a small but steady evolution in the distribution of 
Rb, with all discs becoming progressively larger. Disc scale- 
lengths span today the range 0.8-11 kpc, compared to 0.3-6 
kpc at z ~4. This is due to the inside-out star formation 
scheme adopted in our models (see Sec. 2.1 and Paper II). 
The most probable value has increased only slightly between 
2=4.2 and 2=0, from ~1 kpc to ~1.5 kpc. 

The central surface brightness fio (B) is found to span a 
rather narrow range of values today (20-22.5 mag arcsec -2 ). 
Notice, however, that our models were developed to re- 
produce properties of High Surface Brightness discs, while 
Low Surface Brightness ones (i.e. with fto(B) >22.5 mag 
arcsec -2 ) could be obtained by larger values of the spin pa- 




[Fe/H] [Fe/0] 



Figure 2. Evolution of the distribution functions </>(Q) °f various 
output parameters Q of our models. Results are shown at different 
redshifts, starting with z= 4.3 (dashed curves), then 2=2.8, 1.7, 
1.2, 0.7, 0.5, 0.2 (dotted curves) and z= (solid curves). Top left: 
Distribution of disc scalelength in the B-band Rb ; Top right: Dis- 
tribution of central surface brightness in the B-band Llq.ab(B); 
Second panel left: Distribution of disc restframe absolute magni- 
tude in the B-band Mab(B)\ Second panel right: Distribution of 
disc colour (V-V)ab', Third panel left: Distribution of disc average 
metallicity [O/H]; Third panel right: Distribution of star forma- 
tion rate Bottom left: Distribution of disc average metallicity 
[Fc/H]; Bottom right: Distribution of disc average abundance ra- 
tio [Fe/O]. Our models cover the range 80 < Vc7(km/s) < 360 and 
represent High Surface Brightness discs today (with /j,q(B) < 22.5 
mag arcsec -2 ), but some of them were Low Surface Brightness 
discs in the past. 

rameter A (e.g. Jimenez et al. 1998). We find that, at higher 
redshifts spirals had on average higher central surface bright- 
ness (lower fio(B)), spanning a much broader range of val- 
ues. This is due to the fact that, in our scheme, massive and 
compact discs evolved quite rapidly and thus developed very 
high central surface brightness early on; low mass and ex- 
tended discs evolved slowly, starting from quite low values 
of central surface brightness and developing progressively 
brighter central regions. 

The Ms distribution function was also broader at early 
times, but massive and bright discs became progresssively 
fainter, while low mass ones became increasingly brighter. 
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The most probable Ms value shifted by ~-3 mag in the past 
~10 Gyr. 

The colour distribution function cj>(U — V) shows con- 
siderable evolution in both its shape and average value. As 
time goes on, galaxies become (obviously) redder, while the 
increasing difference in effective ages between massive discs 
(formed early on) and low mass ones (formed much later) 
contributes to broaden the cf>(U — V) function. 

The average metallicity is denned as the amount of oxy- 
gen in the gas divided by the gaseous mass of the disc. It 
increases steadily with time in all discs. At early times, 
massive discs have already developed relatively high O/H 
values, while low mass ones are essentially unevolved; the 
4>(0/H) distribution is then quite broad, as in the case of 
4>(Mb)- At late times, massive discs show little evolution 
while the metallicity of low mass ones increases rapidly, and 
the (j>{0/H) distribution becomes narrower. 

Similar arguments hold for the distribution of the SFR 
4>(^), which spans a large range of values at early times, 
from 10" 3 M yr^ 1 to 60 M yr _1 . At late times, SFR 
values range from 10 _1 M@ yr~ x to 10 Mq yr _1 , with an 
average around 1 Mq yr _1 . 

An interesting consequence of our models is that the 
average Fe/O ratio becomes slightly higher than solar at 
z ~0, since the rate of SNIa (major Fe producers) is such 
that their contribution slighly exceeds the one of O from 
SNII. 

In summary: the distributions of fio(B), Mb, O/H and 
^ become narrower with time and shift to higher average 
values. Both effects result from the rapid late evolution of 
the numerous low mass discs; because of this late evolu- 
tion, the large differences created initially by the rapid early 
evolution of massive discs are reduced at late times. The dis- 
tribution of Rb keeps its overall shape, but its peak value 
slightly increases. Finally, the distribution of (U — V) is the 
only one that broadens with time, because of the increas- 
ing difference between the effective ages of massive and low 
mass discs. 

The results of Fig. 2 are presented in Fig. 3 in a dif- 
ferent form. The evolution of the average values of Rb, 
/tio(-B), Mb , {U — V),0/H and (solid curves) is shown as 
a function of redshift z, along with the ±lcr values of the 
same quantities (inside dark shaded aereas). There is little 
increase in < Rb > for z < 1 (induced mostly by low mass 
discs), while < fj,o > and < $ > decrease slightly at z < 2. 
Other average quantities (< M B >,< U - V >,< O/H >) 
increase constantly, from z ~4 to z—0. In most cases, the 
corresponding values for the Milky Way (dotted curves) are 
larger than the average ones (since the latter are dominated 
by low mass discs). 

In fact, except for the case of Rb and 9, the Milky 
Way values are higher than the average ones by more than 
la during most of the evolution. The Milky Way is a large 
spiral and its evolution is different from the average one 
not only quantitatively, but also qualitatively in some cases 
(e.g. in the case of Mb and fj, where the Milky Way values 
decline for z < 2 while the average ones continue to increase 
or remain constant). 

At this point, we notice that most (but not all) semi- 
analytical models of galaxy evolution today use the Salpeter 
IMF (a power-law with a unique slope -1.35 over the whole 
stellar mass range), which is flatter than the one used here. 




01234 01234 
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Figure 3. Evolution of several quantities computed with our 
models as a function of redshift: scalelength Rb (top left), B- 
magnitude Mab(B) (middle left), average oxygen abundance 
(bottom left), central surface brightness fiQ ^B(B) (top right), 
colour index (U — V)ab (middle right) and star formation rate 
\E> (bottom right). In all panels the solid curve indicates the av- 
erage value of the probability distributions in Fig. 2, while the 
borders of the dark shaded aerea indicate the evolution of the 
corresponding ±lcx values and the light shaded aerea represents 
the remaining discs. The dotted curves indicate the corresponding 
quantities for the Milky Way disc. The dashed curve in the upper 
left figure indicates the predictions of hierarchical models by Mao 
et al. (1998), according to which Rb oc (1 + z)^ 1 . 

In Paper I we discussed at length the reasons of our choice 
and the current observational evidence against a unique 
slope IMF. We also mentioned briefly the main differences 
resulting from its use, mainly the fact that it produces a 
larger effective yield of metals and a smaller number of stars 
of 1-2 Mq, which dominate the galactic luminosity at late 
times. Since the effective yield of a stellar generation de- 
pends not only on the IMF but also on the individual stellar 
yields, which are still uncertain by at least a factor of two 
(see Prantzos 2000 an references therein) , the former effect of 
the choice of the IMF cannot be tested against observations 
today. As for the latter, we find that it results in differences 
of ~30% in luminosity at late times, i.e. ~0.15 magnitudes 
in the B-band; this is also comparable to other theoretical 
uncertainties entering our calculation (neglect of late AGB 
phase in the adopted stellar tracks, crude treatment of dust 
etc.) 



4 EVOLUTION UP TO Z = 1: MODELS VS. 
OBSERVATIONS 

In this section we compare our results to recent observations 
concerning the evolution of discs at moderate redshifts, for 
2 < 1. A nice recent overview of the relevant observations 
can be found in Hammer (1999). 
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Figure 4. Evolution of scalelongth Rb, magnitude Mab(B), 
central surface brightness ^iq^ab(B) and colour index (U — V)ab 
in the range < z < 1.3. The presentation of our results is similar 
to the one of Fig. |i| but the average (solid curves) and ±1<t 
values (borders of the shaded aerea) were calculated taking into 
account only the portion of the distribution function concerning 
large discs (Rb > 47i^" 1 kpc). Observations of large discs by Lilly 
et al. (1998) are shown as circles. Observations of the local sample 
of de Jong (1996) are inside error bars. The Milky Way evolution 
is shown with dotted curves. 



4.1 Large vs. small discs 

Lilly et al. (1998) used two-dimensional surface brightness 
profiles extracted from HST images of galaxies selected from 
the CFRS and LDSS redshift survey to study the evolution 
of several properties of star-forming galaxies between z = 
and z — 1.3. Their observations concern the disc scale- 
length Rb, magnitude Mab(B), central surface brightness 
IJ-ab,o{B) and colour index (U — V)ab- 

The data of Lilly et al. (1998) concern galaxies of var- 
ious inclinations, and are not corrected to the face-on case; 
however, those data do not show any trend with inclination. 
We checked that the scatter produced by varying the incli- 
nation in our models from to 50 degrees is smaller than the 
scatter obtained via the variation of the parameters Vc and 
A. We thus present our results for galaxies seen face-on. The 
comparison to the data of Lilly et al. (1998) with our "face- 
on" models serves merely to check whether the "average" 
history of our models is compatible with the observational 
trends. 

The data of Lilly et al. (1988) concern only large discs, 
with scalelengths Rb >4 h^ kpc. We selected then in our 
models discs satisfying this criterion and we plotted in Fig. 
4 the results concerning the evolution of Rb, fio(B), Mb and 
(U — A)ab in a way similar to the one used in Fig. 3; this 
time, only the part of the distribution functions of Fig. 2 
concerning large discs is used to derive average and ±lcr 
values, and results are shown only in the z <1.3 range, where 
the measurements of Lilly et al. (1998) were made. 

As can be seen in Fig. 4 (left top panel) the data show 
no evolution of i?s up to 2 ~1. The behaviour of our mod- 
els is compatible with this trend, since most large discs are 
already formed at z ~1 and undergo very little increase in 
size afterwards. Note that our models do not reproduce the 
largest galaxies of the Lilly et al. (1998) sample (those with 



Rb > 10 kpc). Such discs are never obtained in our models 
and we consider that they are not representative of "nor- 
mal" spirals like those of the local sample of de Jong (1996) 
on which we based our modelisation in Paper II. We notice 
that the Milky Way has a more rapid late evolution than 
an average disc: its Rb (dotted curve in Fig. 3) decreases by 
a factor of ~2, from ~4 kpc at z=0 to ~2 kpc at z ~1.2. 
In this redshift range the Milky Way behaves almost as ex- 
pected from the hierarchical scenario of Mao et al. (1998), 
i.e. with Rb oc (1 + z) _1 . 

Observed B-magnitudes (left-bottom pannel) exhibit 
some decline since z ~ 1. This is in qualitative agreement 
with our models and is related to the decrease of star forma- 
tion with time in massive galaxies (Paper II). The range of 
Mb predicted by the models is also in fair agreement with 
the observed one. A comparison with Fig. 3 shows that the 
situation differs a lot from the case where all discs are taken 
into account: in that case, the average Mb increases lately, 
because of the late brightening of the numerous small discs. 

There is also a slight decline of the observed central 
surface brightness of discs between z = and z = 1 (right- 
top pannel), in agreement again with the evolution of our 
model cj)(fj,o $ AB(B)) for large discs. 

The right-bottom pannel of Fig. 4 shows the evolution 
of the (U — V)ab colour index. The data seem to indicate 
a strong increase of that colour index with redshift (albeit 
with a large scatter). This is quite at odds with the results 
of our models that predict a slow evolution in this redshift 
range. However, as noticed by Lilly et al. (1998) the colour of 
discs in the local sample of de Jong (1996) has a mean value 
of (U — V)ab ~ 1.5, considerably smaller than their obser- 
vations at moderate redshift (which suggest a mean value of 
2 at z ~ 0.3). Notice that our models reproduce correctly 
colours at redshift z = 0, since they were constructed as to 
match local observations. 

The local values of de Jong (1996), combined to the ob- 
servations of Lilly et al. (1998) point towards little colour 
evolution on average, as in our models. However, the large 
scatter in the (U — V) data can not be explained by our 
histories of spiral galaxies. Brinchmann and Ellis (2000) no- 
tice that colour is a transient property. It could be affected 
by small episodes of interaction, or of enhancement of star 
formation, effects that are not taken into account in our 
models (which give an "average" history of galaxies) and 
could perhaps explain part of the observed scatter. 

The observations of Lilly et al. (1998) suggest a modest 
amount of evolution of large discs in the redshift range < 
z <1, in agreement with the results of our models. This is 
also compatible with the conclusion of Brinchmann and Ellis 
(2000), who found that massive galaxies achieved the bulk 
of their star formation before z ~1. 

Another conclusion of Lilly et al. (1998) is that the main 
evolutionary changes between z = 1 and z = among spirals 
concern small (i.e. with half-light radii (-Re// < 5h^ kpc) 
rather than large discs. Their fig. 17 illustrates the existence 
of small and bright discs at redshift z ~ 1, with no local 
counterparts. Our model also predicts a stronger evolution of 
properties among small galaxies than among large ones, and 
the existence of small and luminous spirals at high redshift. 

This can be seen in Fig. 5, where we present the lumi- 
nosity function 4>(Mb) computed at three redshifts (z=0 on 
the top, z—0.7 in the middle, z=1.2 at the bottom). The first 
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Figure 5. Evolution of the luminosity function of discs as a 
function of size, up to z ~ 1.2. The luminosity function for all 
discs is shown in the first column, while in the second and third 
columns we present the luminosity function of "large" discs and 
"small" discs, respectively. The limit between large and small 
discs is Rb > kpc, adopted from Lilly et al. (1998). In the 

top panels, the functions are shown at redshift ,2=0 (galactic age 
T=13.5 Gyr), in the middle panels at redshift z ~ 0.7 (T=6.5 
Gyr) and in the bottom panels at z ~ 1.2 (T=4.5 Gyr). The local 
(2=0) functions of the top panels (solid curves) are repeated as 
dotted curves in the lower panels, in order to show the evolution 
of (j>(Mg) between low and high redshifts. The local luminosity 
function for all discs of our models (top left panel) is compared 
with two Schechter-typc fits to observational data (grey curves), 
obtained for spiral galaxies with the APM survey (Loveday et al., 
1992) and the SSRS2 survey (Marzke et al., 1998). According to 
our models, large discs were brighter in the past, while small discs 
were both brighter and fainter than today. 



4 
3 

Of 

< 

\ 2 
z 

1 


0.03 



5 



0.01 



OBSERVATIONS 

Local 

..High z 




0.02 





0.03 



Of 
-s. 



0.01 



BRIGHT DISCS 
(M B <-20) 




LARGE DISCS 
(R B >4h 5 - ' kpc) 




-2 -1.5 
= log(R B /V c ) 



-2 -1.5 
Q = log(R B /V c ) 



Figure 6. Evolution of the distribution of Q = log(Rs /Vc), 
where Rg is in kpc and Vc in km/s. Top-left pannel: histograms 
derived by Mao et al. (1998) for a local sample (Courteau 1996, 
1997, solid), and for the high-redshift (< z >~0.7) sample of Vogt 
et al. (1996, 1997; dotted). The comparison of the two samples 
suggests that discs at high redshift were smaller for a given veloc- 
ity. Bottom-left pannel: Distribution computed with our models 
at redshift z=0 (solid curve) and at redshift z ~0.7 (dotted). We 
find a trend similar to the the observationally infered one (albeit 
of slightly smaller magnitude), due to the "inside-out" formation 
of discs in our models. Top-right pannel: Same as bottom-left 
panel, for bright discs only (Mg < —20). Bottom-right pannel: 
Same as bottom-left panel, for large discs only (Rb > kpc). 



brighter by ~2 mag on average in the past. In our models, 
those are moderately massive discs, the external parts of 
which are formed lately. At z ~1 they are already bright, 
but not large enough. This population of "small" and bright 
discs at high redshift has no counterpart today and could 
well constitute the galactic population detected by Lilly et 
al. (1998). We further discuss this point in the end of Sec. 
4.3. 



column presents the results for all discs, the second column 
for "large" discs only (those with Rb > 4/i^q 1 kpc), and the 
third column for "small" discs (Rb < 4/l^" kpc) The local 
luminosity function (all discs, z = 0, top left panel) is com- 
pared with Schechter functions determined for spiral galax- 
ies of the APM survey (Loveday et al., 1992) and the SSRS2 
survey (Marzke et al., 1998). Taking into account that the 
determination of the luminosity function varies from sample 
to sample, that it is affected by the choice of cosmological 
parameters, and that we consider only models of face-on spi- 
ral galaxies, we think that the agreement of the computed 
luminosity function with the observationally derived ones is 
fairly good. 

As already shown in Fig. 4, large discs were on aver- 
age brighter in the past, at least up to z ~1. Fig. 5 reveals 
that they were also less numerous, because some of those 
had Rb < 4/i^" kpc in the past. But the important fea- 
ture of Fig. 5 is that small discs (right panels) were both 
fainter and brighter in the past. A large number of small 
discs (comparable with their observed local population) was 



4.2 The R d - Vc relationship 

The variation of disc size, surface brightness, magnitude etc. 
with redshift can help to test theories of galaxy formation 
and evolution. Mao et al. (1998) used data at low and high 
redshift to deduce that Rd tx (1 + z)~ , in agreement with 
hierarchical models of galaxy formation. 

Before presenting our results and their comparison to 
observations, we notice that the size of a disc is usually 
measured by its exponential scalelength. But a disc can also 
grow by simply increasing its central surface brightness and 
keeping the same scalelength: it is difficult to claim that its 
"size" remains constant in that case (indeed, its R25 radius 
increases) . The use of scalelength to measure size may create 
some problems in that context. 

With this caveat in mind, we proceed to a comparison 
between our models and observation, concerning the quan- 
tity Q=log(i?s /Vc)- The two histograms in the top-left pan- 
nel of Fig. 6 are adopted from Mao et al. (1998). One was 
computed for a local sample of discs (Courteau, 1996, 1997, 
solid), and the other for the high-redshift (z < 1) galaxies 
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Figure 7. Evolution of the distribution of Q = log(Rs) + 
Mg/7.5, where Rg is in kpc. Top-left pannel: histograms derived 
by Mao et al. (1998) for a local sample (Kent, 1985, solid), and 
for the high-redshift (< z >~0.7) sample of Schade et al. (1996, 
dotted). The comparison of the two samples suggests that discs 
at high redshift were smaller for a given magnitude (or, cquiv- 
alently, they were brighter for a given size). Bottom-left pannel: 
Distribution computed with our models at redshift z=0 (solid 
curve) and at redshift z ~0.7 (dotted). We find an effect com- 
patible with the observationally inferred one. However, the ob- 
served high z distribution is broader than ours. Top-right pannel: 
Same as bottom-left panel, for bright discs only (Mg < —20). 
Bottom-right pannel: Same as bottom-left panel, for large discs 
only (R B > 4/1" 1 kpc). 



of Vogt et al. (1996, 1997, dotted). The comparison of the 
two histograms shows that discs in the high redshift sam- 
ple are smaller for a given velocity. In the bottom left panel 
we present our results for all discs. We also find that for a 
given Vc discs were smaller in the past, to an extent slightly 
smaller than the one inferred from observations; in fact, as 
can be seen in Fig. 3, the mean value of our model Rb was 
smaller in the past, up to z ~1, but its rate of variation was 
closer to (1 + z)~ - 5 than to (1 + z)~ x . 

These conclusions could be affected by various bias at 
high redshifts: for instance, the study of Lilly et al. (1998) 
presented in Sec. 3.2 deals only with large and bright discs, 
because fainter discs are not easily resolved. For this rea- 
son, we present in the right part of Fig. 6 the distribution 
of Q=log(7?s/Vc) for bright discs only (Mb < -20, top 
pannel), and for large discs only (Rb > 4Jigo kpc, bottom 
pannel). In both cases, the position of the maximum of <f>(Q) 
is clearly shifted to lower values at high z, but the shift is 
smaller in the case of large discs. Thus, we find that (Rb /Vc) 
should be always smaller on average at high z, but the value 
of the shift depends (albeit not very strongly) on the sample 
selection criteria. 



4.3 The R d — Mb relationship 

Mao et al. (1998) have also investigated the evolution of 
the magnitude-size relationship as a function of z. Fig. 
7 (top-left pannel) shows the distribution of the quantity 
Q=log(i?s) + Ms/7.5 that they obtained for a local sample 



(from Kent 1985) and for the high-redshift data of Schade 
et al. (1996). The latter distribution is shifted with respect 
to the former one towards lower values of Q. From this com- 
parison Mao et al. (1998) deduce that the size- magnitude re- 
lation must have evolved since z ~ 1, in the sense that discs 
of a given luminosity were smaller at early times (or, equiv- 
ently, discs of a given size were brighter). The distributions 
of this same quantity in our models are shown in the left- 
bottom pannel, for redshifts z=0 and z=0.7, respectively. 
We also obtain a shift between the two, compatible with the 
one inferred from the observations. However, our distribu- 
tions (especially the high redshift one) are narrower than the 
observed ones and slightly shifted to larger Q values. The 
reason for this small discrepancy may be the fact that we 
consider discs with Vc >80 km/s, i.e. slightly brighter than 
the lower limit considered in the observed samples. 

The right pannels of Fig. 7 show the same quantity for 
bright discs (top) and large discs (bottom). We find again 
the same effect, i.e. a shift of parameter Q towards higher 
values with time. Selection criteria do not affect then this 
result, namely that discs of a given size were brighter in the 
past. In fact, this is another way of presenting the result of 
Fig. 4, showing that large discs were slightly brighter in the 
past. 

More recently, Simard et al. (1999) analysed 190 field 
galaxies of the DEEP survey in order to determine the 
magnitude-size relationship of discs in several redshift bins, 
up to z ~ 1.1. Their survey is statistically complete for mag- 
nitudes ^814 <23.5. The envelope of their observations at a 
given redshift (taken from their paper) is shown in Fig. ^ 
along with a few other magnitude-size samples (Schade et 
al. 1996, Roche et al. 1998, and Phillips et al. 1997), includ- 
ing the one of Schade et al. (1996) that was discussed in Fig. 
7. The solid diagonal line in each panel is the locus of the 
local Freeman relation hob = 21.65 mag arcsec -2 (Freeman 
1970). To a first approximation, the data follow this relation 
at all redshifts, with some increase of the size dispersion at a 
given magnitude. Also, in the highest redshift bins, the av- 
erage fiQB value of the observed samples is higher than the 
Freeman value. This is another way of presenting the results 
of Fig. 7, namely that at high z discs of a given magnitude 
were smaller, on average. 

Our models of Rb vs Mb are also plotted as a func- 
tion of redshift in the various panels of Fig. 8. We show the 
isodensity contours, as indicated in the bottom-right panel 
of the figure. Our results compare fairly well to the obser- 
vations of Simard et al. (1999) in the redshifts z ~0.2 and 
z -~0.4. At higher redshifts, the dispersion of Rb for a given 
Mb increases, while the average Rb value decreases, both 
in agreement with the observations. In fact, in the highest 
redshift bins, the largest part of our model discs is fainter 
than the cut-off magnitude Mb of the Simard et al. (1999) 
survey; only our brightest discs (the luminous "tip of the 
iceberg") are seen in the observations, and their behaviour 
is entirely compatible with the data. The size dispersion of 
our discs is due to the effect of the spin parameter A on their 
structure and evolution (through the adopted scaling rela- 
tions and the SFR and infall prescriptions, see Sec. 2.2). For 
all A values, massive discs attain their final size early on (be- 
fore z ~1) and evolve little at late times, as discussed in Sec. 
4.1. But, depending on the A value, low mass discs may have 
very different histories: compact ones (low A) evolve rapidly, 
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Figure 8. Disc scalelength Rb vs disc rest frame absolute mag- 
nitude Mg as a function of redshift z. The grey boxes indicate 
various observed samples, with the one of Simard et al. (1999) 
shaded; detection limits are displaced to lower magnitudes as one 
moves to higher redshifts. Our model results, convolved with the 
A and Vc probability functions of Fig. 1, are shown as isodensity 
contours, with the largest number of discs included inside the 
solid curve. At high redshifts, only the brightest of our discs are 
present in the observational data. The grey diagonal line in all 
panels indicates the Freeman value of central surface brightness: 
muB o = 21.65 mag arcsec - 2 (Freeman, 1970). 



i.e. before z ~1, while extended ones (large A) evolve more 
slowly. This explains the large dispersion in the Rb values 
of faint discs at high redshift. 

Simard et al. (1999) draw attention to the fact that the 
apparent increase of mean surface brightness with z in their 
data may be entirely due to selection effects; if the survey 
selection functions are used to correct the surface brightness 
distributions in the different redshift bins, no significant evo- 
lution is found in their data up to z ~1 for discs brighter 
than Mb=-19. 

This is at odds with the results of our models, which 
generically predict an increase of the mean surface bright- 
ness with redshift up to z ~1. This increase is due to the 
fact that most of our discs form their central regions long be- 
fore z ~1; the SFR and the corresponding surface brightness 
of those regions can only decline at late times, i.e. increase 
with redshift. Notice that the situation is different for z >1: 
bright discs exhibit the same trend up to z ~2 (since the 



evolution of their central regions is completed earlier than 
that), while faint galaxies have mean surface brightness de- 
clining with z. According to our senario, they still form their 
inner regions in the z=l-2 redshift range. 

Nine galaxies at z >0.9 of the Simmard et al. (1999) 
sample occupy a region of the Rb — Mb plane rarely pop- 
ulated by local galaxies, i.e. they are bright and relatively 
compact. In our senario, such high surface brightness discs 
do exist at z ~1: they are high Vc and low A discs that have 
completed their evolution earlier than z ~1 and became pro- 
gressively fainter at late times (by ~ 1 mag in Mb in the 
past ~6-7 Gyr. 

4.4 Tully-Fisher relation 

The Tully-Fisher (TF) relation is a strong correlation ex- 
hibited by the whole population of disc galaxies. It relates 
their luminosity L or magnitude M and their circular veloc- 
ity Vc, through L tx y£ /2 ' 5 or M = -a[log{2V c ) - 2.5] + b. 
The exact slope of the relation and its zero-point are still 
the subject of considerable debate (see e.g. Giovanelli et al. 
1997). 

In Paper II we presented a detailed comparison of our 
model results to local (z ~0) TF relations obtained by dif- 
ferent groups in the I and B bands. We noticed that in 
the I band (which presumably reflects better the underlying 
stellar population) the TF slopes of different groups differ 
by more than ~20% (i.e. from a=6.80 in Mathewson et al. 
(1992) to 8.17 in Tully et al. 1998). 

Our results are in much better agreement with the data 
of Han and Mould (the corresponding TF relationship is 
given in Willick et al. 1996) or of Mathewson et al. (1992), 
concerning field spirals, than with those of Tully et al. (1998) 
or Giovanelli et al. (1997), which concern spirals in clusters. 
We notice that our models correspond better to the former 
case than to the latter, since star formation in cluster galax- 
ies may be affected by tidal interactions which are not taken 
into account in our SFR prescription. 

In the B-band, the available local TF relation of Tully 
et al. (1998) concerns again spirals in clusters. Our model 
B— band TF relation is slightly flatter than the observed 
one, but in view of the uncertainties and the caveats men- 
tioned in the previous paragraphs, we do not think that 
this a serious drawback of the models. We note that other 
models (e.g. Mo et al. 1998, Heavens and Jimenez 1999, 
Bouwens and Silk 2000) apparently reproduce better the lo- 
cal B— band TF relation. However, the mass to light ratio 
M/L in those models is either taken as constant or it is 
adjusted as a function of Vc in order to reproduce the TF 
relation. Our own M/L vs. Vc relation (presented in Fig. 8 of 
Paper II) is computed from the fully self-consistent chemo- 
spectrophotometric evolution models, which reproduce all 
the main properties of local discs without any adjustment 
except for the initial scalings (see Sec. 2.2 and Paper II). 

Using spectroscopic observations of the Keck telescope 
and high-resolution images from the Hubble deep field, Vogt 
et al. (1996, 1997) determined a Tully-Fisher (TF) relation 
for 16 galaxies lying between redshifts z=0.15 and 1 in the 
rest-frame B— band. They find no obvious change in the 
shape or slope of the TF relation with respect to the local 
one. Assuming then the same slope, they derive a modest 
brigthning of 0.36T0.13 mag between redshifts z ~0 and 
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Figure 9. Evolution of the Tully-Fisher relation as a function of 
redshift z (from top to bottom) in the B— and I— bands (left and 
right, respectively). The average relationship obtained with our 
models is shown as solid curve, among ±lcr values (dark shaded 
aerea) and all remaining discs (light shaded aerea); in the bottom 
three panels, the dotted curve is the average model TF relation 
of the top panel (z=0), allowing to measure the evolution. Boxes 
in all the left panels represent the observed local TF relation 
in the B— band (Tully et al. 1998), while boxes in all the right 
panels show the local TF relation in the 7— band (Mathewson et 
al. 1996). Data at high redshift in the left panels are from Vogt 
et al. (1996, 1997). 



1. Their data, split in two redshift bins (0.2< z <0.6 and 
0.6< z < 1, respectively) are shown in Fig. 9. 

We also present in Fig. 9 the TF relationship computed 
with our grid of models at four galactic ages, corresponding 
to z ~0, 0.4, 0.8 and 1.2, respectively (from top to bottom). 
In the top panel, results are compared to the local B— band 
observations of Tully et al. (1998). In the second and third 
panels, our model TF relations are compared to the data of 
Vogt et al. (1996, 1997). Our results are certainly compat- 
ible with the data, but in view of the large error bars and 
selection effects (only the brightest galaxies are detected in 
the the z ~0.8 bin) no strong conclusions can be drawn. It 
is clear, though, that our models predict a steepening of the 
TF relation in the past, since massive discs were brighter, 
on average, than today, while low mass discs were fainter. 
This steepening with repect to the local slope is small and 
difficult to detect at z ~0.4 and more important in higher 
redshifts. We also notice that at high redshifts the disper- 



sion in luminosity becomes larger for the fainter discs. In 
both cases, obervations of faint galaxies (Mb >-18) would 
be required to establish the behaviour of the TF relation at 
high redshift. 



5 EVOLUTION AT REDSHIFT Z >1. 
5.1 Scalelength distributions 

The size-luminosity relation at redshifts z >1 has been re- 
cently investigated by Giallongo et al. (2000). They used 
data from the Hubble Deep Field-North, where morpholog- 
ical information is available for all galaxies up to apparent 
magnitude / ~26. Using colour estimated redshifts and as- 
suming exponential scalelengths, they derived absolute blue 
magnitudes Mb and disc scalengths Rd as a function of red- 
shift. They compared these data, concerning the redshift 
range 1< z <3.5, to those of the ESO-NTT deep field, 
concerning the redshift range 0.4< z <1 (Fontana et al. 
1999). The results of their study appear in Fig. 10 (his- 
tograms) for "bright" and "faint" discs (upper and lower 
panels, respectively). Comparing to predictions of "stan- 
dard" semi-analytical models, including merging histories 
for dark haloes, Giallongo et al. (2000) found a significant 
excess (by a factor of ~3) in both redshift ranges of bright 
and small discs (i.e. discs with Rd <2 kpc in the low red- 
shift range and with Rd <1.5 kpc in the high redshift range); 
they also found a smaller excess (by ~50%) of "faint" discs 
in both redshift ranges. 

Our results are also plotted in Fig. 10. In each panel, 
dotted and dashed curves indicate scalelength distributions 
obtained at the two extremes of the corresponding redshift 
range. When adding the model distributions of the upper 
and lower panels at a given z one obtains the total distribu- 
tion of Rb for that z, shown in Fig. 2. We notice that the 
observed and model distributions are normalised differently 
(corresponding scales are indicated on the left and right axis 
of each panel) , but the comparison of their shapes allows to 
draw some interesting conclusions. 

In the low redshift range (left panels), our model distri- 
butions evolve very little between 2=1. and z=0.4, with both 
"bright" and "faint" discs becoming larger with time. In 
both cases, distributions at z=l compare slightly better to 
the data than distributions at z—0.4. The overall agreement 
with the data is excellent. We find no defficiency of small 
discs (either bright or faint ones), contrary to Giallongo et 
al. (2000). We notice that we always use the probability dis- 
tributions Fy and Fx of Fig. 1, assuming no merging. This 
assumption turns out to be crucial in the agreement of mod- 
els with the data. Indeed, taken at face value, the results on 
the left panel of Fig. 10 imply that the adopted prescriptions 
for the scaling relations, infall, SFR, Fv and Fx are sufficient 
to make discs with the correct distributions of scalelengths 
at intermediate redshifts. 

In the high redshift range (right panels of Fig. 10) we 
also find little evolution and excellent agreement with ob- 
servations for the "faint" discs (Mb >-20). We notice that, 
by construction, our models predict little early evolution for 
discs of low mass and luminosity. For "bright" discs, we find 
significant evolution in the Rb distribution, in the sense that 
this distribution is significantly shifted to larger and more 
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Figure 10. Distribution functions of the scalelength in the B— band Rb in two different redshift ranges: 0.4< z <1. {left) and 
l.< z <3.5 (right), respectively. In both redshift ranges, distributions concern "bright" discs (upper panels) and "faint" discs (lower 
panels), respectively; notice that the magnitude limits are different in the two cases (Mg=-19 on the left and Mg=-20 on the right). 
Observed distributions (histograms) are from the surveys of the ESO-NTT deep field at low z and of the Hubble Deep Field at high z 
(Giallongo et al. 2000); the corresponding scales on the vertical axis are displayed on the right of each panel. In each panel we show two 
model distributions (in dotted and dashed curves, respectively) corresponding to the two extremes of the corresponding redshift range, 
i.e. at 2=0.4 and 1 in the left panels and at z=l and 3.5 in the right panels, respectively. Our model distributions are not normalised to 
the observed ones and the corresponding scales are displayed on the left of the vertical axes. 



numerous discs when the redshift decreases from z=3.5 to 
z=l: in that redshift range a large number of initially small 
and "faint" discs (belonging to the lower right panel of Fig. 
10 at 2=3.5) become brighter and larger at z=l (and "mi- 
grate" to the upper panel), due to the adopted inside-out 
formation scheme. However, none of the two model distri- 
butions compares well to the observed histogram, and this 
will obviously be true for any other model distribution at 
redshift 1< z <3.5: we find then (as Giallongo et al. 2000) 
a serious deficiency (factor ~3) of small and bright discs 
at high redshift. We notice, though, that this is the only 
discrepancy of our models with respect to the observations, 
while Giallongo et al. (2000) found discrepancies in all four 
panels of Fig. 10. 

The defficiency of our models concerns bright and small 
discs, i.e. discs with high surface brightness. Selection effects 
may play some role in that case, since they obviously favour 
the detection of such discs, at the expense of more extended 



ones, having lower surface brightness. If there is no selection 
bias in the data, then it is the first hint that our models 
(based on the simple prescriptions of Sec. 2) may have some 
difficulties in reproducing observables of the high redshift 
Universe. In fact, our models are "calibrated" on the Milky 
Way disc and the infall timescale adjusted as to reproduce 
all major observables of discs in the local Universe (Papers 
II, III and IV). The dicrepancy encountered in Fig. 10 may 
imply that mergers between small discs (not taken into ac- 
count here) did play some role in that redshift range: such 
mergers would affect only a small fraction of the population 
of the lower right panel of Fig. 10, but would contribute sig- 
nificantly in populating the deficient part of the curves in 
the upper rigt panel. 
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Figure 11. Upper panels: Relation between the abundance gra- 
dient of oxygen dlog(0/H)/dR and the inverse of the scalelength 
Rg as a function of redshift. The relationship found for local discs 
in Paper III (Prantzos and Boissier 2000) is supported by obser- 
vations. Data at z ~0 from Garnett et al. 1997 (open squares) 
and van Zee et al. 1998 (filled circles with error bars); the one 
dicrepant point of van Zee et al. (1998) concerns an extremely 
small and low surface brightness disc. We find a similar anticor- 
relation also at higher redshifts; discs evolve "homologuously" . 
Lower panel: The abundance gradient of our models, expressed 
in dex/iig and weighted by the probability functions Fy and F\ 
of Fig. 2, is plotted as a function of redshift (Solid curve: average 
value; dark shaded aerea: ±lcr values; light shaded area: all other 
discs. 



5.2 Abundance gradients vs. scalelengths 

In Paper III of this series (Prantzos and Boissier 2000) we 
showed that a remarquable correlation exists between the 
abundance gradient of a disc (expressed in dex/kpc) and its 
scalelength Rb- small discs have large (in absolute value) 
abundance gradients and vice versa. This is easy to under- 
stand qualitatively, since abundance gradients are created by 
radially varying quantities, like the SFR txR -1 adopted 
here. In a small disc (say, Rb=2 kpc), such quantities vary 
significantly between regions which are distant only 1 kpc 
from each other; a large abundance gradient is obtained in 
that case. In a large disc (say Rb~6 kpc) regions separated 
by 1 kpc have quite similar SFR; such discs cannot develop 
important abundance gradients, in dex/kpc. 

Although intuitively obvious, this property has never 
been properly emphasised before. In Paper III it was shown 
that there is a unique relation between dlog(0/H)/dR and 
1/ Rb for all discs. This relation is found in observed nearby 
spirals and is fairly well reproduced by our models of discs at 
2 ~0, as can be seen on the upper left panel of Fig. 11. The 
importance of that correlation is obvious: the abundance 
gradient can be derived if the disc scalelength is known and 
vice versa. 



In this section we explore whether this correlation is 
valid at higher redshifts. In Fig. 11 we plot the relation of 
dlog(0/H)/dR vs. 1/Rb for our model discs at redshifts z~0 
(where they compare quite favourably with observations), 
z=1.5 and z~3. At all redshifts we find the relation obtained 
locally. The "homologuous evolution" of spirals, invoked by 
Garnett et al. (1997) for local discs, is also found in the high 
redshift Universe. 

In the lower panel of Fig. 11 we present then the quan- 
tity dlog(0/H)/di?s of our models (i.e. the abundance gra- 
dient expressed in dex/Rs, instead of dex/kpc). We plot this 
quantity, weighted by the probability functions Fv and F\, 
as a function of redshift. It can be seen that during most of 
the history of the Universe, this quantity remains ^constant, 
at a value of dlog(0/H)/dRs ~-0.20 dex/Rs; only at 
very late times, this quantity increases slightly (in absolute 
value), been dominated by the numerous small discs that 
evolve considerably in the z ~0.-0.5 range. For that same 
reason, the dispersion around this value increases at low z. 
Both the stability of the relation of dlog(0/H)/dR vs. 1/ Rb 
and the decrease in the dispersion of dlog(0/H)/di?s with 
redshift are novel and important predictions of our model. If 
true, they will allow to infer something about the chemistry 
of high redshift discs from their morphology. 

At this point, we notice that Damped Luman-a sys- 
tems (DLAs) may be (proto-)galactic discs such as those 
considered here: in a previous paper (Prantzos and Boisser 
2000) we have shown that, when observational bias are taken 
into account, the observations of Zn abundance in such sys- 
tems as a function of their redshift (0.5< z <3.5) are nicely 
explained by our models. The apparent "no chemical evolu- 
tion" picture suggested by the observations of those systems 
(e.g. Pettini et al. 1999) is quite compatible with our under- 
standing of the chemical evolution of discs. 



6 SUMMARY 

In this paper we explore the implications of our disc galaxy 
models for the high redshift Universe. Our models (presented 
in detail in Papers I, II, III and IV) are "calibrated" on the 
Milky Way, utilise simple prescriptions for the radial varia- 
tion of the SFR and infall rate and simple scaling relations 
(involving rotational velocity Vc and spin parameter A) and 
reproduce all major observables concerning discs in the local 
Universe. A crucial ingredient for the success of our models 
is the assumption that massive discs form the bulk of their 
stars earlier than their low mass counterparts. 

We assume here that the probability functions Fy and 
Fx of the two main parameters of our models are indepen- 
dent and do not evolve in time, i.e. we assume that merg- 
ers never play a significant role in the overall evolution of 
disc galaxies. We follow then the evolution of the distri- 
bution functions of various quantities and we compare the 
results with available observations concerning discs at high 
redshifts. Our results are summarised as follows: 

1) Most large discs (Rb >4hJ 1 kpc) have basically com- 
pleted their evolution already by z ~1. Subsequently, they 
evolve at ~constant scalelength, while their luminosity and 
central surface brightness decline (Sec. 4.1). These features 
of our models are in quantitative agreement with observa- 
tions of the CFH survey (Lilly et al. 1998) . Our models also 
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predict a mild evolution of U— V colour of large discs in that 
redshift range, while the Lilly et al. (1998) data suggest a 
slightly stronger evolution; but their data, extrapolated to 
z ~0, are in disagreement with the local values derived by 
de Jong (1996), while our models reproduce fairly well the 
observations of the local sample. Since U — V colours can 
be strongly affected by "transient" phenomena (e.g. short 
"bursts" of star formation), we do not think that this dis- 
crepancy is a real problem for the models. 

2) Because of the "inside-out" formation scheme of discs 
adopted in our models, we predict that discs were, on av- 
erage, smaller in the past for a given rotational velocity Vc 
or for a given magnitude Mb- Both results are in qualita- 
tive agreement with recent samples of discs in redshifts up 
to z ~1 (Sec. 4.2 and 4.3). Our models reproduce fairly well 
the data of Simmard et al. (1999) on the Rb vs. Mb relation 
as a function of redshift, up to z ~1. In particular, our mod- 
els produce some compact and bright discs that Simmard et 
al. (1999) find at z ~1 and which have no counterparts in 
the local Universe (Se. 4.3). 

3) We find that the Tully- Fisher relation was steeper 
in the past, with massive discs been, on average, brighter, 
and low mass discs been fainter than today. Our models are 
in agreement with observations of the TF relation in the 
B— band obtained by Vogt et al. (1997) at intermediate 
redshifts (z ~0.4 and 0.8), taking into account their error 
bars (Sec. 4.4). However, these data are insufficient to probe 
the evolution of the TF relation, which requires data on 
fainter discs (down to Mb ~-18 in that redshift range). 

4) Our models reproduce failry well the distribution of 
disc scalelengths of the ESO-NTT survey, concerning both 
bright and faint discs in the 0.4-1 redshift range (Sec. 5.1). 
They also reproduce well the scalelength distribution of faint 
discs in the 1-3.5 redshift range obtained in the Hubble Deep 
Field (Giallongo et al. 2000). On the contrary, we fail to 
reproduce the correpsonding distribution of bright discs in 
that same redshift range: our models show a deficiency of 
small (Rb ~1 kpc) and bright (Mb <-20) discs with re- 
spect to the Hubble Deep Field data. If this discrepancy is 
not explained in terms of selection effects (favouring the de- 
tection of those high surface brightness discs) it could be a 
strong indication that merging of small discs played indeed 
an important role in that redshift range. 

5) The anticorrelation between the abundance gradient 
dlog(0/H)/dR (expressed in dex/kpc) and the scalelength 
Rb, that we found in Paper III for local discs, is found to 
be valid also at higher redshifts (Sec. 5.2). We find that the 
abundance gradient of discs, expressed in dex/_R_B, varies 
very little with redshift and presents a smaller dispersion 
at high redshifts. These findings establish a powerful link 
between the optical morphology and the chemical proper- 
ties of disc galaxies, valid in both the local and the distant 
Universe. 

In summary, we explored the consequences of our "back- 
wards" model of disc galaxy evolution for the high redshift 
Universe. We found no serious discrepancies with currently 
available data up to redshift 2=1. At higher redshifts, we 
find hints that the model may require some modifications, 
provided that currently available data are not affected by 
selection biases. 
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